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ABSTRACT 

We examine the a-element abundance ratio, [a/Fe], of 5620 stars, observed by the Sloan Extension 
for Galactic Understanding and Exploration survey in the region 6 kpc < i? < 16 kpc, 0.15 kpc 
< \Z\ < 1.5 kpc, as a function of Galactocentric radius R and distance from the Galactic plane \Z\. 
Our results show that the high-a thick disk population has a short scale length (Lthick ~ 1-8 kpc) 
compared to the low-a population, which is typically associated with the thin disk. We find that the 
fraction of high-a stars in the inner disk increases at large \Z\, and that high-a stars lag in rotation 
compared to low-a stars. In contrast, the fraction of high-a stars in the outer disk is low at all \Z\, 
and high- and low-a stars have similar rotational velocities up to 1.5 kpc from the plane. We interpret 
these results to indicate that different processes were responsible for the high-a populations in the 
inner and outer disk. The high-a population in the inner disk has a short scale length and large scale 
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height, consistent with a scenario in which the thick disk forms during an early gas-rich accretion 
phase. Stars far from the plane in the outer disk may have reached their current locations through 
heating by minor mergers. The lack of high-a stars at large R and \Z\ also places strict constraints 
on the strength of radial migration via transient spiral structure. 



1. INTRODUCTION 

Detailed observations of the Galactic thick disk can 
be used to constrain the relative importance of cos- 
mological accretion and secular processes in the forma- 
tion and growth of the Milky Way d i sk. T hick disk 
stars are old (> 8 Gyr, iBensbv et al.l I2005D and pro- 
vide a fossil record of the Galaxy at z ^ 2. Obser- 
vations of these old stars can serve as a complement 
to studies of distant galaxies at early times, many of 
which are seen to be thick, turbulent, cl umpy, star- 
forming disks (e.g., Elmcgr een fc Elmegre cn 2005, 120061 : 
iForster Schreiber et al.. ,20091 |2011[ ). Observations of 
nearby galaxies have show n that thick disks are common , 
with similar properties (IDalcanton &: BernsteinI 120021 : 
lYoachim fc DalcantonI [20051 l2006l l2008bllaD . which sug- 
gests that thick disks are a generic feature of disk galax- 
ies. Thus, the processes responsible for the existence of 
the Milky Way thick disk may play an important role in 
the formation of all disk galaxies. 

Several mechanisms have been proposed to ex- 
plain the formation of the thick disk. Four such 
mechanisms are: (1 ) vertical hea t ing through mi- 
nor mergers (e.g., K azantzidis et al.l 120081: iRead et al] 
2008t iVillalobos fc Hclmi .2008; iKazantzidis et al.l 120091: 



Purceil et al.l 120091: iBird et al.H2012D: (2) dire ct accre- 
tion of stars from satellites (jAbadi et al.ll2003D : (3) in- 
situ formation during an early tur bulent disk phase 
due t o high gas accretion rat es (e.g., iBrook et al.l 12004 
I2005t iBournaud eFal] I2009D : and (4) radial migra- 
tion of stellar orbits via reso nant interactions wi th 
transient spiral stru cture ( e.g.. iRoskar et al] l2008bl lal: 
ISchdnrich fc Binnev 2009a, bl: ILoebman et al.ll201lD . Sce- 
narios 1-3 fit within the context of hierarchical structure 
formation as predicted by ACDM cosmology, while sce- 
nario 4 is possible in a disk in complete isolation. Each of 
these scenarios can be tested through comparisons with 
the observed chemical and kinematic properties of stars 
in the Milky Way. 

Sin c e the thick disk's disco very by star counts (jYoshiil 
119821 : IGilmore fc ReidI I1983D . an apparent dichotomy 
between thin and thick disk populations has been 
established. Stars belonging to the thick disk are 
older and more met a l poor (e.g. IGilmore et al.l [l995l: 
Chiba fc BeeTi [2OO0I: IBensbv et al.l 12004 llvezic et al.l 
20081) . In addition, thick disk stars have chemical 



abundance patterns distinct from the thin disk, with 
thick dis k stars being enhanced in a- and r-process el 
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The results of these latter studies show that the scatter 
in the observed [a/Fc] is small, suggesting that thick disk 
stars formed quickly in a well-mixed interstellar medium. 

The high [a/Fe] characteristic of thick disk stars indi- 
cates that they formed in a period of rapid star formation 

yyc@ucolick.org 
^ Packard Fellow 



and chemical enrichment, during which Type II SNe were 
able to contribute significant amounts of a-clcmcnts into 
the interstellar medium before Type la SNe increased the 
abundance of iron-peak elements. The observed abun- 
dance trends of nearby thick disk stars have been used to 
estima te that the th ick disk formed o ver a peri od of ^1- 
3 Gvr (iGratton et al. 2000: Mash onkina & Gehren 200fl 
iMashonkina et al.ll2003l: IBensbv etral.ll2004[) . The chem- 
ical properties of thick disk stars are thus a powerful tool 
for understanding the chemical enrichment and star for- 
mation history of the Galaxy. 

Whether the enhanced [a/Fe] for thick disk stars in 
the solar annulus (i?GC.0 = 8.0 kpc) is present for m- 
situ thick disk stars (i.e., those at large \Z\, where the 
thick disk is expected to dominate) remains an open 
question. Ea rly in-situ studies ex amined only the [Fe/H] 
distribution (jGilmore et al.l[T995| ). Analyses of more el- 
ements has largely been restricted to stars in the so- 
lar neighborhood, where stars are typically divided into 
thin and thick disk popu lations by t heir kinematics (e.g. , 
IBensbv et ai]|200ll2005l ). Recentlv. IBensbv et al.l (|201lf ) 
examined the abundances of a sample of 119 red giants 
in the inner and outer disks using high-resolution spec- 
troscopy and found evidence that the [a/Fe] trends ob- 
served in the inner disk do no t extend to the outer dis k 
(see also IBensbv et all l2010al: lAlves-Brito et al.l I2010D . 
Thus, while many of the thick disk fo rmation scenarios 
disc ussed above (such as the wo rk of iBrook et al.ll2005l 
and iSchonrich fc BinnevI l2009b( ) are able to reproduce 
the dichotomy seen in the chemical properties of thin 
and thick disk stars in the solar neighborhood, observa- 
tions have only begun to test the models at a wide range 
of R and |.^| in the Galactic disk. 

Using a sample of old disk stars from the Sloan 
Extension for Galactic Unde rstanding and Exploration 
(SEGUE; lYannv et al.ll200l survey, we have begun to 
explore a larger volume of the Galaxy. We previously 
showed that the radial metallicity gradient in [Fe/H] be- 
comes flat for stars located at vertical heights \Z\ > 1 
kpc from the Galactic plane, where the thick disk is 
expec ted to be the dom inant population (jCheng et al.l 
I2OI2I hereafter iPaper J ) . This result is consistent with 
a chemically homogeneous thick disk, which is predicted 
by thick disk formation during a period of early gas-rich 
accretion (scenario 3). The flat gradient could also be 
explained if the strength of radial mixing is sufficient to 
e rase a pre -existing gradient in the thin disk (scenario 4). 

IPaper J also demonstrated that the reported flatten- 
ing trend in the radia l metallicity grad ient at i? > 10 
kpc (|Yong et al.l 120051: iLuck et al.ll2006D could arise be- 
cause all of the distant tracers were located at large \Z\, 
where the radial gradient is flat. We found that the ob- 
served trend could result from a simple superposition of 
a negative radial gradient at small \Z\ with a flat ra- 
dial gradient at large \Z\. Because the clusters discussed 
in the literature at large R were also located at large 
|Z|, it is unclear whether the observed trends are due to 
changes in the radial or vertical directions. Therefore, 
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we stressed the importance of examining the trends in 
[Fe/H] as a function of both R and \Z\. 

In this work, we extend our analysis of abundance gra- 
dients in the Milky Way to examine the a-element abun- 
dance ratio, [a/Fe], as a function of both Galactocentric 
radius R and distance from the plane \Z\, using a sam- 
ple of 5620 field sta rs from SEGUE. Our work is comple- 
mentary to that of iBensbv et alj ()2011|) . as our sample 
of main sequence turnoff stars is more than an order of 
magnitude larger than their samp le of 119 K giants. We 
use a subset of the sample from iPaper J . which covers 
the region 6 kpc < i? < 16 kpc, 0.15 kpc < \Z\ < 1.5 
kpc. We present our data and results in [21 and [21 re- 
spectively. In []4l we present estimates for the thin and 
thick disk scale lengths; the procedure and errors are de- 
scribed in more detail in the Appendix. We discuss the 
implications of our results in ij5l 

2. DATA 

Our sample consists of main sequence turnoff stars 
from the Sloan Extension for G alactic Unde r stand- 
ing and Exploration (SEGUE, lYannv et all 120091: 
lEisenstein" et al. 2 0TT1). part of the Sloan Digital Sky Sur- 
vey (SDSS; lYork et al.l l 2000[). The data are o btained us- 
ing the same CCD camera (jGunn et al.lll998D. telescope 
(iGunn et al.l |2006[) , and filter system (jFukugita et al.l 
11996D as the SDSS. In this paper, we use a subset of stars 
from th e sample of 7605 main sequence turnoff stars in 
IPaper J . which cover the region 6 kpc < i? < 16 kpc, 0.15 
kpc < \Z\ < 1.5 kpc. Briefly, these stars are selected us- 
ing a cut in g — r color. Stellar parameters Tog, log 
[Fe/H], and [a/Fe] are determined from low resolution 
(i? - 2000) spec tra using the SE G UE Stellar Parameter 
Pipel ine (SSPP.'Lee et al."2 008allbl: lAllende Prieto et al.l 
I2OOI .Smolinski et al. 2011i: ILee et al.l I2011b[ )! In the 
present work, we select the stars with sufficient signal-to- 
noise (S/N > 20 pixel""'^, where each pixel corresponds 
to ^ lA) for good [a/Fe] measurements, which yields 
5771 stars with [a /Fe] measured to a precision of 0.1 dex 
(ILee et al.ll20TTbh . 

The S/N cut effectively imposes a magnitude limit, 
and because the [a/Fe] of a star will not affect its mag- 
nitude significantly, this magnitude limit does not bias 
our sample against low-a stars. Dartmouth isochrones 
(jDotter et al.l 120081 ) show that a 0.2 dex difference in 



[a/Fe] is equivalent to < 0.1 mag difference in g-band 
magnitude, which corresponds to a 10% error in distance. 
This is a small effect compare d to the total error in dis- 
tance 20% — 25% estimated in IPaper^ . Therefore, we do 
not expect any significant systematic differences between 
the volumes sampled by high- and low-a stars. In our 
sample, the high- and low-a stars span the same ranges 
in R and \Z\, and they are seen out to the same distances. 

As described in IPaper^ . we assign a weight to each tar- 
get to correct for selection effects. The weight accounts 
for three properties of the selection: (1) Objects in re- 
gions with the highest extinction in each line of sight were 
not considered for spectroscopy. (2) Not all candidates 
for spectroscopy are observed. (3) The g — r color cut 
introduc es a bia s against redder, metal-rich stars. We 
show in IPaper II that using the calculated weights suc- 
cessfully reproduces the true metallicity gradients in a 
mock catalog (§6.3.2). For a detailed discussion of the 
selection biases and how we correct for them using our 



weighting scheme, see §4 and the Appendix of IPaper J . 

In addition to the weights we calculated in IPaper J . we 
apply a weight to account for the stars that are removed 
by the additional S/N cut imposed on this sample; this is 
a small effect compared to the other weights. Taking all 
stars with non-zero weights, we are left with a sample of 
5620 main sequence turnoff stars. Most of the stars that 
are given weights of zero are very blue objects, which are 
likely to be hotter stars that are not on th e main se quence 
(for more discussion see the Appendix of iPaperT ). 

Distances were determined using photometric paral- 
lax methods, by comparing the SSPP stellar p arameters 
and ph otometry to the theoretical isochrones of lAn et"al] 
()2009D . as described in IPaper 1 . Usi ng a mock catalog of 
stars g enerated from the model of iSchonrich &: BinnevI 
(j2009a| ). we estimate t he errors in the distances to be 
~ 20% - 25%; see §6 of IPaper J for details. 

3. RESULTS 

3.1. Abundance Trends as a Function of R and \Z\ 

Figure 111 shows [a/Fe] as a function of Galactocentric 
radius R, in four slices of \Z\, for our sample of main 
sequence turnoff stars, color coded by [Fe/H]. Most of 
the high-a population is confined t o small radii {R < 1 
kpc), consistent with the results of IBensbv et al.l ()201lD . 
who found a lack of high-a stars in the outer disk. In 
our sample, this lack of high-a stars is seen at all \Z\. 

Figure 121 shows abundance trends in [Fe/H] and [a/Fe]. 
T he top left panel shows the solar neighborhood sample 
of IBensbv et all (|2003ll2005l) . in which stars are assigned 
to the thin and thick disks (red and blue, respectively) 
according to their kinematics. In the solar neighbor- 
hood, kinematically hot stars (i.e., thick disk stars) are 
a-enhanced relative to kinematically cooler stars (i.e., 
thin disk stars) at the same [Fe/H]. The top right panel 
shows the total SEGUE sample, where we see two pop- 
ulations analogous to the solar neighborhood thin and 
thick disk stars: (1) low-a stars that, like solar neigh- 
borhood thin disk stars, appear to follow a linear trend, 
with [a/Fe] slightly decreasing as [Fe/H] increases, and 
(2) a tail of high-a stars that, like solar neighborhood 
thick disk stars, are more metal poor than the low-a 
population. 

In the remainder of our analysis, we divide our sample 
into high- and low-a stars, with the goal of comparing the 
high-a (low-a) stars to the kinematically selected thick 
(thin) disk stars in the Bensby sam ple; this i s simila r 
to the chemical separation done by ILee et ahj ()2011aD . 
We make the cut at [a/Fe]= -1-0.2, where the number of 
stars appears to drop dramatically at large i?, as seen in 
Figure 111 This is marked by the horizontal dotted line 
in the top right panel of Figure |2l The number of high- 
and low-a stars is indicated in the top right corner in 
parentheses. 

Finally, in the bottom four panels, we show the SEGUE 
sample divided into four R,\Z\ bins with cuts at i? = 10 
and \Z\ = 0.5 kpc. The numbers in the top right corner 
of each panel indicate the weighted fractions of the high- 
and low-a populations, with the raw number of stars in 
parentheses. The weighting has the effect of slightly de- 
creasing the fraction of high-a stars in each bin, but the 
effect is not dramatic. This is likely because sele ction on 
g — r color is biased against metal-rich stars (see IPaper D 
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Fig. 1. — a-element abundance ratio [a/Fe] vs. Galactocentric radius R in four \Z\ slices. The SEGUE data are shown as dots, colored 
coded by [Fe/H]. At all \Z\, the majority of the high-a stars are located at small R (< 10 kpc). 



for details); metal-poor stars, which are more likely to 
be high-a stars, are weighted less heavily to compensate 
for the bias. 

At small R (< 10 kpc, left panels), we see the same 
high- and low-a populations as in the total SEGUE sam- 
ple. The presence of two populations is especially evident 
at \Z\ > 0.5 kpc. The fraction of high-a stars increases 
toward large \Z\ (> 0.5 kpc), from 12% to 31%, which is 
what we expect if the contribution from a high-a thick 
disk is greater far from the plane. At large R (> 10 
kpc, right panels), in contrast to what is seen at small 
R, the fraction of high-a stars docs not increase at large 
\Z\ (> 0.5 kpc); the fraction is low at all \Z\. This ob- 
servation at large R is inconsistent with the picture of 
a high-a population associated with a thick disk compo- 
nent, which should become more dominant at large \Z\ 
at all R. The main result of Figure [2] is that the [Fc/H]- 
[a/Fe] properties for stars at small R are consistent with 
those found for solar neighborhood stars, while at large 
R there is a lack of high-a stars. Furthermore, the frac- 
tion of high-a stars at large R does not increase with \ Z\, 
as expected if there is a high-a, thick disk population in 
the outer disk. 

The change in fraction of high-a stars at large R sug- 
gests that the chemical abundances of stars change with 
Galactocentric radius, even at large \Z\. But how can 
this result be reconciled with the flat rad ial meta llicity 
gradient at \Z\ > 1.0 kpc that we found in iPaper^ Fig- 
urc[3]shows the radial metallicity gradient A[Fe/H]/Ai? 
in four \Z\ slices for the total sample (gray, left column) 



and divided into low-a (pink, middle column) and high- 
a stars (blue, right column). For each sample, we fit a 
linear trend to the data, with e ach point weighted to ac- 
count for selection biases, as infPaQeO- The number of 
stars and the slope of a linear fit to the data arc indicated 
in the bottom right corner of each panel. 

The change in A[Fe/H]/Ai? with \Z\, as shown in Fig- 
ure [Sj is summarized in Figured The radial gradient of 
the high-a sample (blue triangles) is flat at all R and \Z\, 
but it is not solely responsible for the flattening trend 
with \Z\ seen in the total sample (gray squares). The 
high-a stars do, however, make the gradients for the total 
sample flatter, especially at \Z\ > 1.0 kpc, where the frac- 
tion of high-a stars is the largest. The flattening trend of 
the low-a sample (pink diamonds) is closely followed by 
the trend in the total sample. The results of lPaper^ are 
also shown (black circles) in Figure ID These are slightly 
different than the gradients measured for the total sam- 
ple because of the S/N cut imposed on the sample in 
this work, but are still within the uncertainties. The er- 
rors are estimated using 500 Monte Carlo realizations of 
the data, where we perturb the stellar parameters Toff, 
[Fe/H], and [a/Fe] by the typical er rors (20 K, 0.3 dex 
and 0.1 dex, respectively); see §6 of iPaper^ for details. 

3.2. Kinematics of the High- and Low-a Populations 

In addition to different chemical properties, thin and 
thick disk stars in the solar neighborhood exhibit dif- 
ferent kinematic properties. A comparison of the kine- 
matics of high-a stars at large and small R can help 
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Fig. 2. — Abundance trends [a/Fe] vs. [Fe/H]. Top left panel: the solar neighborhood sample of lBensbv et al.l II2003I . |2005| ). with thin 
and thick disk stars (red and blue, respectively) assigned according to their kinematics. Top right panel: the total SEGUE sample. The 
horizontal dotted line indicates where we make the distinction between low- and high-a stars. Bottom four panels: the SEGUE sample, 
divided into four bins of R and \Z\. The labels on the contours indicate the number of objects in a box with dimensions of 0.15 dex in 
[Fe/H] and 0.05 dex in [a/Fe]. In each panel, the weighted fraction of high- and low-a stars (blue and red, respectively) is indicated, with 
the raw number of stars in each population in parentheses. The abundance patterns at < 10 kpc arc similar to those seen for thin and 
thick disk stars in the solar neighborhood, with the fraction of high-a stars increasing at large \Z\. At ij > 10 kpc, the fraction of high-« 
stars is low at all \Z\. 



distinguish whether high-a stars at large R are the outer 
disk tail of the inner disk population, the high-a tail 
of the outer disk population, or a different population 
altogether. In this section we examine the rotational ve- 
locities Vcf, of high- and low-a stars as a function of R 
and \Z\. For this analysis, we only consider 3985 stars, 
which have good proper motions, as described below. 

We calculate three-dimensional velocities in Carte- 
sian coordinates {U,V,W), and polar coordinates 
{Vb.,V^,Vz), using radial velocities along with proper 
motions o btained by combin ing the USNO-B and SDSS 
catalogs (IMunn et al.l I2004D . We use the criteria of 
IKilic et all (|2006l) " to obtain a sample with clean proper 
motions: sigRa < 525 mas, sigDec < 525 mas, 

match = 1, nPit = 6, dist22 > 7", where sigRa 
and sigDec are the residual for the proper motion fit in 
Right Ascension and Declination, match is the number of 
objects within a 1" radius, nPit is the number of plates 
the object was detected on, and dist22 is the distance 
to the nearest neighbor with g < 22. The efficacy o f 
these criteria have been explored bv lDong et al.l (|201lD . 

These selection criteria identify 3985 stars in our sam- 
ple with reliable proper motions; this subsample has the 
same distributions in distance, magnitude, color, [Fe/H], 
and [a/Fe] as the larger sample of 5620 stars, so we treat 



it as a representative sample and do not apply any addi- 
tional weights. Statistical errors on the proper motions 
are roughly 3 — 3.5 mas yr""'^, which corresponds to a 
tangential velocity error of 28 — 33 km s^^ at a distance 
of 2 kpc. 

In this section we present histograms, corrected using 
our weighting scheme, which are needed to properly ac- 
count for the different sampling along different lines of 
sight. For example, there are many more stars in the 
lines of sight toward smaller R {I < 90°) compared those 
at the anti-center {I ~ 180°), but the number of spectra 
are approximately equal in all directions. Accounting for 
this effect is necessary to reproduce the correct distribu- 
tions, but doing so magnifies the Poisson noise. 

Figure [5] shows histograms for the Bensby sam- 
ple (top left panel), the total SEGUE sample (top right 
panel), and the SEGUE sample in the four R,\Z\ bins 
(bottom four panels), similar to Figure [2l In each panel, 
as in Figure [21 the sample is divided into high- and low- 
a stars (blue and red, respectively) at [Q;/Fe]= 0.2. The 
shaded regions indicate errors estimated by generating 
500 Monte Carlo realizations of our data, where we per- 
turb the stellar parameters, radial velocities, and proper 
motions; the typical errors on the radial velocities and 
proper motions are 6.0 km s~^ and 3 mas yr~^, respec- 



6 



1 


-1 
-2 

1 


-1 
-2 

1 


-1 
-2 

1 


-1 
-2 



N = 364 
Slope: 0.00045 



N = 251 
Slope: -0.024 



1.0 < IZI < 1.5 



N = 113 
Slope: 0.0093 



X 




0.50 < IZI < 1.0 



N = 2333 
Slope: -0.017 



N = 1806 
Slope: -0.045 



N = 527 
Slope: 0.0045 



X 





N = 2212 
Slope: -0.065 



N = 1856 
Slope: -0.082 



0.25 < IZI < 0.50 



N = 356 
Slope: -0.015 



X 



0.15 < IZI < 0.25 



Total 
N - 550 

Slope: -0.11 



[a/Fe] < 0.2 
N - 469 

Slope: -0.14 



^^^^^^^ 



[a/Fe] > 0.2 
N - 81 

Slope: -0.0070 



10 15 
R (kpc) 



205 



10 15 
R (kpc) 



205 



10 15 
R (kpc) 



20 



i 



Fig. 3. — Galactocentric radius R vs. [Fc/H] in four \Z\ slices for high- and low-o stars. Left column: Total sample. Middle column: 
Low-o stars ([a/Fe] < +0.2). Right column: High-a stars {[«/Fe] > +0.2). In each panel the raw number of stars and the measured slope 
are indicated in the bottom right corner. The lines show a linear fit to the data, with each star weighted to account for selection biases. 
The spacing of the symbols on the linear relation indicates the radial distribution of the stars. The radial gradient in the total and low-a 
samples become flatter at large \Z\, while the radial metallicity gradient of the high-a sample is flat at all \Z\. 

each population (vertical solid lines) using an outlier- 
resistant algorithm, which trims values greater than 
three median absolute deviations from the median (ver- 
tical dotted lines in Figure [5]) . The numerical values of 
(V^) are indicated in the top right corner of each panel, 
with the raw number of stars in parentheses. The errors 
on the means are indicated by the shaded regions sur- 
rounding the vertical solid lines. For the low-a samples, 
the error on the mean is 2 — 3 km s~^, while for the 
high-a samples, which have fewer stars, it is ~ 4 — 10 km 
s-i. 

At i? < 10 kpc, the kinematics of high-a stars are con- 
sistent with those seen for thick disk stars in the solar 
neighborhood. Figure [5] shows that at all \Z\, they lag in 
rotation behind low-a stars by ~15 km s""'^, in rough 
agreement with measurements in the li t erature (e.g., 
Chiba fc Beersl[2000t ISoubiran et al.l[200l ICarollo et ahl 
2010[) . which have values ~ 20 — 50 km s~^, depending on 
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Fig. 4. — Radial metallicity gradient, A[Fe/H]/AiJ, vs. distance 
from the plane, \Z\, for high- and low-o stars. The radial gradients 
for the t otal sam ple in this work (gray squares) are consistent with 
those ofl^jgeFJ, which were measured using a larger sample with 
a less strict S/N cut (black circles). The change in the radial 
gradient of the total sample as a function of \Z\ is driven by the 
change in the radial gradient of the low-o stars (pink diamonds); 
the radial gradient of the high-a stars (blue triangles) shows no 
obvious trend and is consistent with a flat gradient (slope of zero) 
at all \Z\. The number of stars used in each gradient measurement 
in the low- and high-a samples is indicated. The error bars reflect 
the random errors in the gradi ent meas urement due to errors in 
the stellar parameters (see §6 of IPapeFTI for details). 



tively. 

We calculate the mean rotational velocities (V^) 



for 



how the thin and thick disk populations are separated. 
Like thick disk stars in the solar neighborhood, the high- 
a stars at i? < 10 kpc also have wider distributions in Vz 
and Vr. This result implies that they belong to a kine- 
matically "hotter" population, which has larger random 
motions in the radial and vertical directions in addition 
to a larger lag in V^. For reference, the mean, median, 
median absolute deviation, and skewness of the popula- 
tions' distributions in V^, Vz, and Vr for both high- and 
low-a stars arc tabulated in Table [TJ The errors on (V^) 
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R < 10 
IZI > 0.5 



196.6 (628) 
180.5 (321) 



R < 10 
IZI < 0.5 



202.9 (1392) 
187.9 (269) 




[o/FeJ < 0.2 (3212) 
[a/Fe] > 0.2 (773) 



R > 10 
IZI > 0.5 



218.5 (854) 
210.8 (129) 
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Fig. 5. — Weighted rotational velocity distributions for high- and low-a stars. Top left panel: the solar neighborhood sample of 
IBensbv et al.l II20031 . 120051 ). with thin and thick disk stars (red and blue, respectively) assigned according to their kinematics. Top right 
panel: the SEGUE sample, with high- and low-o stars (blue and red, respectively) divided at [o/Fe]= -1-0.2. The shaded regions indicate 
the errors estimated using 500 Monte Carlo realizations of our data. Bottom four panels: the SEGUE sample, divided into four bins of 
R and \Z\. In the top right corner of each panel, the mean rotational velocity (V^) of each population is indicated, with the raw number 
of stars in each population in parentheses. To calculate an outlier- resistant (V^) (vertical solid lines), we exclude targets that are more 
than three median absolute deviations from the median value (vertical dotted lines). The errors on (V^) are indicated by the surrounding 
shaded regions. At < 10 kpc, the high-a stars lag in rotation behind the low-a stars by ~15 km s~^. At R > 10 kpc, the difference in 
(V^) between high- and low-o stars is < 8 km s~^, which is within the measurement errors. The different kinematic properties of high-a 
stars at large and small R suggest that they may be different populations with different origins. 



are indicated. 

In addition, there is a large fraction of high-a stars 
with low V0 (< 150 km s~^); this tail is not present 
in the low-a population. One way to quantify this fea- 
ture is to examine the skcwness of the distributions 
of both populations. We find that for the high-a popu- 
lation, the skcwness of the distribution is larger than 
for the low-a population; at \Z\ < 0.5 kpc {\Z\ > 0.5 
kpc) the distribution of the high-a population has a 

skcwness of -1.77to;2 (-1.27lo 3), "^^ilc ^0 distribu- 



4-0.3 



tion of the low-a population has a skcwness of — 0.50_q 

(-o.2ot°:3)- 

The skewed shape of the distribution of both high- 
and low-a stars at i? < 10 kpc is consistent with a pop- 
ulation that is falling in density with increasing R. If 
stars that have slow (fast) rotational velocities are on or- 
bits with guiding centers within (beyond) the solar circle, 
then they are interlopers fr om the inner (outer) disk (e.g., 
IBinnev fc Merrifieldlll998[ ) ; the skcwness of the distribu- 
tion results from the higher stellar densities in the inner 
disk compared to the outer disk. Thus the larger skew- 
ness of the high-a population in is consistent with 
there being a steeper density gradient for high-a stars 



compared to that of low-a stars. The steeper density 
gradient of the high-a stars is consistent with the pic- 
ture of the high-a population having a short scale length 
compared to the low-a population. 

At i? > 10 kpc, however, high- and low-a stars have 
the same mean rotational velocities, within the errors, 
and the fraction of high-a stars with low is compa- 
rable to that of the low-a stars. In addition, the widths 
of the Vz and Vn distributions of high- and low-a stars 
are similar. While at i? < 10 kpc, the different distri- 
butions are indicative of two populations with different 
structural parameters, at i? > 10 kpc, we are unable to 
distinguish between the kinematics of high- and low-a 
stars. If high- and low-a stars at i? > 10 kpc are part of 
the same population, then the observations imply that 
high-a stars at large and small R may have different ori- 
gins. 

4. SCALE LENGTH OF THE HIGH-« POPULATION 

The lack of high-a stars at large R suggests that the 
high-a population, which is typically associated with the 
thick disk, has a short radial extent. The similarity be- 
tween the V,f, distributions of high- and low-a stars at 
large R also supports this idea. In this section, we es- 
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TABLE 1 

Properties of Velocity Distributions 



R < 10, \Z\ > 0.5 




R > 10. \Z 


> 0.5 














Vz 




Va 








Vz 




Vr 




Low- a 


High-a 


Low- a 


High-a 


Low- a 


High-Q 


Low- a 


High-Q 


Low-Q 


High-a 


Low-Q 


High- Of 


Mean 
Median'' 


196.6+^;^ 
198.8 


180.5+-*;;, 
181.5 


-2.4 
-0.6 


-6.0 
-8.1 


108.2 
106.0 


93.9 
93.3 


218.5+^'^ 
217.8 


210.8+»3^ 
208.7 


18.8 
16.6 


21.8 
20.7 


-34.0 
-30.5 


-37.2 
-40.4 


mad'''" 


25.2 


35.3 


24.5 


41.5 


37.0 


43.8 


29.6 


34.8 


33.3 


31.0 


31.9 


35.5 


Skew 




_1 77+0-3 

'-■ ' ' -0.2 


0.08 


0.34 


0.20 


-0.42 






-0.25 


-0.07 


1.77 


0.23 


R < 10, \Z\ < 0.5 




R > 10. \Z 


< 0.5 














Vz 




Va 








Vz 




Vr 




Low-Q 


Higli-a 


Low- a 


High-Q 


Low- a 


High-Q 


Low- a 


High-Q 


Low-Q 


High-a 


Low-Q 


High-a 


Mean"'" 
Median'' 


202.9+J 3 
203.0 


187.9l*j 
192.0 


1.1 
2.3 


-3.9 
-6.5 


50.0 
52.9 


46.1 
46.3 


226.3l-^« 
224.9 


233. 4tJ^;^ 
240.7 


16.1 
15.1 


18.2 
11.7 


-41.1 
-38.6 


-25.6 
-10.4 


MAD''-'^ 


19.2 


34.1 


17.0 


24.9 


35.8 


45.2 


18.7 


20.3 


20.1 


15.9 


22.6 


23.2 


Skew 


-0.20t°j 


-1 27+"* 

J^--^' -0.3 


-0.13 


-0.04 


-0.18 


-0.80 


-0.33+Si 




0.06 


-0.07 


-0.28 


-0.98 



Outliers clipped at three median absolute deviations 
'' Mean, median, and MAD reported in km s~'. 
" Median Absolute Deviation 



timate the scale length of the a-enhanced thick disk by 
quantifying the fraction of low- and high-a stars as a 
function of R and \Z\. We then compare the data to 
expected values based on different combinations of thin 
and thick disk scale lengths, ithin and ithick- In addi- 
tion, we require that the predicted total stellar density as 
a function of R and |.^| is consistent with th e total stel- 
lar de nsity of the best-fit model obtained by lJuric et aD 
(|2008D . While it is possible that the data may be better 
fit by different analytical functions, we restrict our anal- 
ysis to a radial exponential profile for the disk, for which 
we also have total stellar density measurements. 

In using the fractions of high- and low-a stars to esti- 
mate the thin and thick disk scale lengths, we implicitly 
assume that high-a stars belong to the thick disk and 
low-a stars belong to the thin disk, as is observed in the 
solar neighborhood. We include the second constraint 
that t he total stellar density match that of I Juric et alJ 
(|2008l l because our data provide a poor constraint on 
the local normalization of the disk models we fit. In our 
fits, w e fix the sca l e heig hts and the total normalization 
to the I Juric et"al] ()2008D values. We exclude the lowest 
\Z\ slice, where our sample does not cover a large range 
in R. In this section, we present the results of our analy- 
sis; details about the our procedure and the uncertainties 
in our estimates are provided in the Appendix. 

The blue shading in Figure [5] represents the reduced 
values we obtain by comparing the observed and ex- 
pected values of the high- and low-a fractions, as a func- 
tion of R and | Z | , for different combinations of ithin and 
-^thick- The thin yellow contours show the discrepancy 
between the total stellar density predicted by each com- 
bination of scale lengths and the total densit y from the 
best-fit scale lengths reported by I Juric et "all ()2008l see 
also their Figure 21). Com bining both the co nstraints 
of total stellar density from I Juric et "all (120081 ) and the 
high- and low-a fractions from our sample, we find the 
best combination of thin and thick disk scale lengths to 

be Lthin = 3.4^0.9 kpc, Lthick = l-Stoi kpc, marked by 
the large orange cross in Figure [6] The best fit values 
and error bars are obtained by marginalizing over each 
axis and determining the 68% confidence interval for each 
scale length. The 68% contour in two dimensions for the 
combined constraint is shown as the thick dashed orange 
line. More details about the procedure and error analysis 



are given in the Appendix. The blue-shaded map shows 
that for any given thin disk scale length, the preferred 
thick disk scale length, as constrained by the high- and 
low-a fractions (i.e., the white regions), is always shorter 
thin disk scale length. 

We note that our measurement of the thin and 
thick disk scale lengths is not well constrained, as we 
do not have data at large radius near the midplane. 
Our preferred value for the thin disk scale length is 
slightly larger, but consistent with values based on near- 
infrare d photometry fro m the Spacelab Infrared Tele- 
scope (jKent et al.l Il991 | . ^ 3.0 kpc) and the Cosmic 
Background Explorer (jFreudenreichI 119981 ~ 2.6 kpc). 
Other SDSS/SEGUE studies, which use detailed model- 
ing to estimate disk structural parameters, are also con- 
sistent with these values for the t hin disk scale length 
(jJuric et al.ll2008l: iBovv et al.ll20lH) . While our sample's 
spatial coverage is not ideal for constraining the thick 
disk scale length, it is an improvement on earlier stud- 
ies. Furthermore, our best-fit value d oes not change sig - 
nificantly when we vary the assumed I Juric et "all ()2008[ ) 
scale lengths or remove possible halo contaminants. See 
the Appendix for more details. 

Figure [7] shows the weighted fractions of high- and low- 
a stars (blue and red squares, respectively), as a function 
of R in three \Z\ slices. The shaded regions indicate the 
errors on the measured fractions. Each column shows 
the same data compared to the expected fractions of thin 
(pink dotted lines) and thick (blue solid lines) disk stars 
for three different combinations of thin and thick disk 
scale len gths: (1) L^u = 2-6, ithick = 3.6 kpc in the left 
column (jJuric et alJ |2008D (2) Lthin = 3 8, Lt niek = 2.0 
kpc in the middle column (iBensbv et al.ll20l"ll) . and (3) 
-^^thin = 3.4, -Lthi ck =1-8 kpc i n the right column (this 
work). For the iBensbv et al.l (1201 If ) scale lengths, we 
show the predictions for both a constant thick disk scale 
height (thin, dark blue line) and one that varies as Rg / L 
(thick, light blue line), as described in their paper. 

As seen in Figures [1] and [21 Figure [7] shows that the 
fraction of high-a stars decreases at large R, in every \Z\ 
slice. If we assume that the thick dis k is populat e d onl 
by high-a stars, then the model of I Juric et "ahl (|200' 



> 



vastly overpredicts the fraction of high-a stars a,t \Z 
0.5 kpc; in our bin i? > 10 kpc, \Z\ < 0.5 kpc (FigureE]) 
we would expect ^ 50% of our sample 700 stars) to be 
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Fig. 6. — Best-fit thin and thick disk scale lengths. The blue shading shows the reduced values using the fractions of high- and low-a 
stars as a function of R and \Z\ as a constraint. The thin yellow contours show the constraint of the total stellar density as determined 
by lJuric et al.l II2008I) . Our best estimate of the thin and thick disk scale lengths using both constraints are Ltjjin = 3.4]to'g kpc and 
^thick = 1-8^0 5 '^P'^' rnarked by the large orange cross. The 68% contour for the combined constraint is shown as the thick dashed orange 
line. The published values of lJuric et all ||2008|) and IBensbv et all l l20Tll'l arc indicated in green and purple, respectively. 



enhanced in [a/Fe] instead of the ^ 10% (~ 200 stars) 
that we observe. The data agree much better with a 
sh orter thick disk scale length, consistent with the results 
of IBensbv et al.l (|2011l . middle column). Using [a/Fe] as 
a proxy for membership in the thin and thick disks in 
our larger sample, we estimate that the thick disk has a 
shorter scale length than the thin disk (right column). 

5. DISCUSSION 

5.1. The Thick Disk Scale Length 

The results presented above show that the fraction of 
high-a stars drops off at large R and that the high-a 
populations at small and large R have different kinematic 
properties. Both of these results a re consistent with th e 
properties of the K-giant sample of IBensbv et al.l ()20lH ). 
who found that the lack of high-a stars was consistent 
with the thick disk having a shorter scale length than 
the thin disk (Lthin = 3.8, ithick = 2.0 kpc). Using our 
data, we estimate the scale lengths to be Lthin = 3.4lQ'g 

and Lthick = l-8!to 5 '^P'-- While the scale lengths are 
not well constrained with our data, the thick disk scale 
length is consistently found to be shorter than 2 kpc (see 
the Appendix). 

Independent analyses with other SDSS/ SEGUE sam- 
ples have also found simil ar scale lengths for stars with 
thick disk chemistry. iCarollo et al.l (|2010( ). using the ve- 



locity ellipsoid of stars in a particular metallicity range 
and location in the disk (-0.8 < [Fe/H] < -0.6, 1 < 
\Z\ < 2 kpc), estimated the thick disk sca le length to be 
2.2 kpc. In recent work. lBovv et al.l ()201lD measured the 
scale length of their "a-old" population (—1.5 < [Fe/H] 
< -0.25,0.25 < [a/Fe] < 0.50) to be 1.96 kpc. While 
the methods differ, these studies reach the same conclu- 
sion: the population of stars associated with a thick disk 
component in the solar neighborhood has a short radial 
scale length. 

It is worth emphasizing, however, that the above re- 
sults apply for stars with particular properties. The 
results presented in this work reflect the fractions 
of high- and low-a stars, and our scale length es- 
timate reflects the radial extent of the high-a pop- 
ulation, which we associate with the thick disk 
based on studies of the solar neighborhood. Previ- 
ous s t udies of exter nal galaxies ([Dalcanton fc Bernstein! 
I2002t lYoachim fc Dalcanton ,2001) and the Milkv Wav 
(|Juric et al.l 120081: Ide Jong et al.M2010t) have relied on 
surface brightness profiles and star counts, respectively, 
which follow the total stellar density, with no information 
about the stellar populations. 

We now introduce terminology to distinguish be- 
tween these two methods of identifying the two com- 
ponents of the disk. First, we will refer to the struc- 
tural thin and thick disks to describe the components 
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Fig. 7. — Weighted fractions of high- and low-a stars (blue and red squares, respectively) vs. Galactocentric radius R in three \Z\ slices. 
The shaded regions indicate the errors on the measured fractions. Pink dotted and blue solid lines show the expected contributions of 
the thin and thick disks, respectively, for different combinations of scale l engths: (1) -Lthin = 2.6, ithick = 3.6 kpc in the left column 
llJuric e t al. 2008), (2) 

^thin = 3.8, Lthick = 2.0 kpc in the middle column I IBensbv et al.ll2011l) . and (3) Lthin = 3.4, Lthick = 1-8 kpc in 
the right coluimi (this work). In the middle column, we also show the expected fractions if the scale height varies linearly with Rq/L, as 
discussed by Bcnsbv ct al. (2011, thinner lines). Our data support a shorter scale length for the high-Q population, in agreement with the 
conclusion of .Bensbv et al. (.20 11 ). 

the scale lengths of structural thick disks seen in exter- 
nal galaxies. Observationally, structural thin and thick 
disk scale lengths h ave been found to be uncorr elatcd 
in external galaxies (jPalcanton fc Bernstein|[2002f ). For 



that are identified using the total stellar densities, ei- 
ther thro ugh star counts or s urface brightness pro- 
files (e.g., iGilmore fc ReidlflQSSl : iDalcanton fc BernsteinI 
[200l) . Second, we will refer to the chemical thin and 
thick disks to describe the components that are identi- 
fied using the chemical and/or kinematic properties o f 
stars (e.g., iBensbv et all [200l [2005t ILee et al.l IMTal ). 
In our work, we have found that the chemical thick disk 
has a shorter scale length than the chemical thin disk. 

Whether star counts in the SEGUE imaging along our 
lines of sight support a short scale length for the struc- 
tural thick disk remains an open question. Previous es- 
timates of the scale length of the structural thick disk 
in the Milky Way have typically relied on star counts 
in higher latitude data and therefore do n ot have sig- 
nificant leverage in the r adial direction (e.g.. lJuric et alj 
120081: iChang et alll20lH ). Our lines of sight reach larger 
R at small | Z \ and may provide additional constraints on 
the scale lengths of the structural thin and thick disks. 
Comparing star counts in our low latitude lines of sight to 
different combinations of structural parameters, as well 
as exploring different radial profiles, will be the focus of 
future work. 

If the structural thick disk is found to have a short 
scale length, in agreement with our result for the chem- 
ical thick disk, then we can compare these results to 



galaxies with circular velocities greater than ^100 km 
s~^, structural thin disks with larger scale lengths than 
structural thick disks h ave b een reported (see Table 5 of 
lYoachim &: Dalcanton, 20061 ). Having a structural thick 
disk with a short scale length, then, would not make the 
Milky Way an unusual galaxy. 

5.2. Possible Relation to the Hercules Thick Disk Cloud 

In Table [2l we list the weighted fraction of high-a 
stars along each of the 11 lines of sight in our sample. 
Most of the 11 lines of sight have fewer than 15% of 
their stars with [a/Fe] > -1-0.2. Three lines of sight at 
i? < 10 kpc, however, have ~ 20 — 50% of their stars with 
[a/Fe] > -1-0.2. These lines of sight are directed toward 
the Hercules Thick Disk Cloud, a stellar o verdensity in 
the di s k, wh ich h as been studied in detail bv lLarsen et al.l 
(l2OT0l [2OTTI) andlHumphrevs et al. ( 2011.) . 

iLarsen et al 1 ([201(1 )' detect this overdensity as an excess 
in star counts in the first quadrant (Ql, 0° < I < 90°), 
compared to the fourth quadrant (Q4, 270° < / < 
360°), in particular, at Galactic coordinates 20° < / < 
55°, 20° < 6 < 45°. Stars associated with this over- 
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TABLE 2 

Fraction of High-« Stars per Line of Sight 
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0.305 


2554 


2564 


94.0 


14.0 


600 


165 
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2672 


187.0 
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81 


0.121 


0.094 


2678 


2696 


187.0 


8.0 


580 


67 


0.116 


0.080 


2712 


2727 


203.0 


8.0 


546 


93 


0.170 


0.142 



density in Ql lag in rotation by 30 km compared to 
stars in Q4, but ha ve metallieities simi l ar to stars in anal- 
ogou s fields in Q4 (jParker et alJl200l : [Humphreys et alJ 
[20Tl[) . Their preferred scenario for the existence of the 
overdensity is that dynamical interactions with the bar 
cause stars to pile up in a "gravitational wake" (e.g., 
iHernauist fc Weinberdll992D . This feature may be re- 
lated to the Hercules-Aquila Cloud see n in the SDSS 
(|Belokurov et al.ll2007t lJuric et al.ll2008| ). 

The rotation rates, w, in our three high-a lines of sight 
(22 — 31 km kpc~^) are slightly larger than the Ql 
fields of lHumphrevs et all ()201lL 15 - 26 km s'^ kpc'^), 
which is not unexpected, since our lines of sight are at 
lower Galactic latitude. A direct comparison is not pos- 
sible because our samples do not overlap spatially. Ad- 
ditionally, our sample does not cover a sufficiently large 
part of the Galaxy to fully test the presence of an asym- 
metry: all of our inner disk stars are in Ql, and we have 
no stars in Q4 with which to make a comparison. We 
do not currently have the necessary data to confirm or 
exclude the possibility that these three lines of sight are 
associated with this overdensity. 

5.3. Implications for Thick Disk Formation 

The short scale length of the chemical thick disk can 
be used to constrain various scenarios for thick disk for- 
mation, such as the four described in SJT] In the following 
discussion, we will use the generic term "thick disk" to 
refer to both the structural and chemical thick disks. We 
assume these to be the same, as is done in all of the mod- 
els discussed. 

The lack of high-a stars at i? > 10, \Z\ > 0.5 kpc 
puts an upper limit on the strength of migration due 
to transient spiral structure (scen ario 4) ; the iV-body 
simulation of lLoebman et all ()2011t) , for example, predict 
that this mechanism can transport many high-a stars 
from the inner disk to large R and \ Z\. In this simulation, 
high-a stars are present at all R because they are old 
and have had more time to migrate to large R. The lack 
of high-a stars that we have observed at large R, then, 
implies that this mechanism cannot be very efficient. 

If the stars we observe at large R and \Z\ reached 
their current positions through radial migration, the lim- 
ited extent of the high-a population could be evidence 
that the mechanism must have some radial dependence 
on the strength of migration. One such mechanism re- 
lies on the presence of a bar and a steady state spiral 
pattern, which leads to more efficient mixing at cer- 
tain radii (e.g.. iMinchev fc Famaevll2010l : iBruiietti et al.l 



120 lit ). iBrunetti et"an (|2011| ). for example, find that stars 
that are close to the corotation radius of the bar are 
more likely to migrate. If the high-a stars in our sam- 
ple were born in the bulge and have since migrated to 
where we observe them at i? < 10 kpc, this would explain 
the similar abundance patterns t hat have been reported 
for thick disk and bulge stars (iMel endez et al.f 120081: 
Bensby et al.ll2009|: lAlves-Brito et alE oiO: Be nsbv et atl 
2010bHRyde'eraLll2010l:IGonzalez et al.i,2011i ). Oneway 
to test this scenario is to examine whether the kinematic 
properties of stars in these simulations are different for 
those mixed in the inner and outer disks. 

If the thick disk does indeed have a shorter scale length 
than the th in disk, as is suggest ed by both our data and 
the data of IBensbv et al.l (j2011|) . and radial migration is 
not the dominant mechanism, our results may have impli- 
cations on the formation and merger history of the Milky 
Way disk. A range of thick dis k scale lengths can result 
from different merger histories. iBrook et al.l ([20p . [2001 
showed that chaotic gas accretion at early times (sce- 
nario 3) results in a thick disk with a shorter scale length 
than the thin disk, while an early gas- rich 2:1 merger 
of two disks results in a thick disk with a longer scale 
length. The variation in merger histories would provide 
a possible explanation for the range of structural thin 
and thick disk scale lengths observed in nearby galax- 
ies ([Dalcanton fc BernsteinI 120021 : lYoachim fc DalcantonI 

\mm . ^ 

Scenarios involving minor mergers can also be con- 
strained. If thick disk stars originated from an initially 
thin disk (scenario 1) then any heating event must have 
occurred in a primordial disk with a short scale length. 
The predominantly low-a stars at large R and \Z\, then, 
should come from a more chemically evolved disk. Radial 
mixing induced by late minor mergers has b een shown 
to be very efficient for stars in the outer disk (iBird et al.l 
I2012D and could explain the presence of low-a stars at 
large R and \Z\. However, predictions can be depen- 
dent on the particular models being examined and their 
initial conditions (e.g.. IdT Matteo et al.l[20Tl[ ). If thick 
disk stars originated from an accreted satellite (scenario 
2), then stars contributed by a single satellite should be 
located in a torus-like structure, and we should see the 
same abundance trends in Ql and Q4. Our current sam- 
ple is insufficient in spatial coverage to explore this pos- 
sibility. 

Finally, the radial gradients in [Fc/H] for high- and 
low-a stars (sec Figures [3] and [3]) also provide constraints 
on thick disk formation. Radial migration due to tran- 
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sient spiral structure (scenario 4) could explain the flat- 
tening trend in the low-a stars, but the mechanism is 
too efficient in current simulations, resulting in too many 
high-a stars in the outer disk, especially at large \Z\. The 
observed distributions of [Fe/H] and [a/Fe] could be ex- 
plained by the following: First, early gas-rich mergers 
(scenario 3) created a chemically homogeneous, high-a 
population in a thick disk with a short scale length. Sub- 
sequently, a low-a thin disk forms and is heated by minor 
merger activity at later times (sc enario 1), mixing stars 
at large R (e.g., iBird et ani2012[ ) and flattening the ra- 
dial metallicity gradient at large \Z\. 

6. SUMMARY 

We have demonstrated, using a sample of 5620 main 
sequence turnoff stars from the SEGUE survey, that the 
high-a population, which is associated with the thick 
disk in the solar neighborhood, has a short scale length 
(ithick ~ 1.8 kpc) and a flat metallicity gradient at all 
\Z\. The abundance trends at i? < 10 kpc show a di- 
chotomy between high- and low-a stars similar to that 
seen between thick and thin disk stars obse rved in the 
solar neig hborhood (jBensbv et al.]l2003L l2005f ) . The frac- 
tion of high-a stars increases with \Z\, and these high-a 
stars lag in rotation compared to low-a stars (by ~15 
km s~^), similar to the difference in kinematics seen 
for thin an d thick disk star s in the solar n eighborhood 
(|Chiba fc B eers 2000,; Sou biran et al.ll2003| ). 

At i? > 10 kpc, the fraction of high-a stars is lower 
than at small R and does not increase with \Z\. High-a 
stars at large R also do not lag in rotation compared to 
low-a stars, with both populations having similar mean 
rotational velocities. These results suggest that the high- 
a stars in the outer disks may simply be the tail of the 
[a/Fe] distribution; the stars far from the plane in the 
outer disk {R > 10, \Z\ > 0.5 kpc) likely have different 
origins than those far from the plane in the inner disk 
(i? < 10, \Z\ > 0.5 kpc). 

The fractions of high- and low-a stars are consistent 
with the expected values fo r a thick disk with a short 
scale length as suggested by iBensbv et al.l ()2011|) . Using 
the fractions of high- and low-a stars as a function of 
R and \Z\, we estimate the thick disk scale length to be 
ithick 1.8 kpc. In addition, it is possible that the lines 
of sight in our sample with large fractions of high-a stars 
are related to the Her cules Thick Disk Cloud, a stellar 
overdensity studied by iHumphrevs et al.l ()201lD . A sam- 
ple of stars with better spatial coverage, particularly in 
Q4, is required to fully explore this possibility. 

We find that the presence of a thick disk with a short 



scale l e ngth is consistent with the scenario of lBrook et al.l 
()2004 I2005D . in which the thick disk formed during a 
turbulent disk phase at early times when gas accretion 
rates were high. In the outer disk, stars may have been 
moved to large R and \Z\ through radial mixing due to 
late minor mergers (e.g., iBird et aI1l2012[ ). The lack of 
high-a stars can be used to constrain the strength of 
radial migration of stars from t he inner disk induced 
by transient spiral structur e ( e.g.. iRoskar et ah 2008b. 3; 
iSchonrich fc BinnevI l2009airbl : iLoebman et"anT2011ll . If 
stars in the outer disk arrived at their current locations 
through radial migration, some radia lly-dependent mech- 
anisms may be responsib le (e.g., iMinchev fc FamaevI 
IMotlBT^etti et al.llMl . 
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APPENDIX 
APPENDIX: SCALE LENGTH ESTIMATES 
PROCEDURE 

To determine the scale lengths of the thin and thick disks, we follow the prescription of lJuric et al.l (|2008| ) and model 
the Galaxy as two double exponential disks plus a two-axial power-law ellipsoid halo (their Equations 21-24): 

p{R, Z) ~ Pd{R, Z; ithin, -f^^thin) 

+/pD(i?, Z; Lthick, i^thick) + Ph{R, Z) (B1) 

where 

Rr. 



pM Z; L, H) ^ Pj,{Rq,0) X e^p [ 
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TABLE 3 ^_ 

Structural Parameters Measured by IJuric et al.I POOSI ) 



Parameter Bias-Corrected Error Definition 



Zq 


25 pc 


20% 


Solar offset from the Galactic plane 




2600 pc 


20% 


Thin disk scale length 




300 pc 


20% 


Thin disk scale height 


S 


0.12 


10% 


Thick disk normalization relative to thin disk at i? — Rq , Z — 


-t'thick 


3600 pc 


20% 


Thick disk scale length 


-f^thick 


900 pc 


20% 


Thick disk scale height 


/h 


0.0051 


25% 


Halo normalization relative to thin disk at R — Rq , Z — 


gn 


0.64 


;S 0.1 


Halo ellipticity 


riH 


2.77 


< 0.2 


Halo power law 



R Z + Zq 



and 



(B2) 



PH(i?,^) =PD(i?0,O)/H 



Rc 



(B3) 



The definitions of the parameters and the values used (the bias-corrected parameters in their Table 10) are listed in 
Table [3l We vary the thin and thic k disk scale l ength s and fix the remaining parameters (including the scale heights 
and the total normalization) to the lJuric et al.I (|2008[ ) values, which are constrained using their photometric sample. 
Our spectroscopic sample is not well suited for determining the total stellar density because of the smaller sample size 
and the pencil-beam nature of the observations. 

For each combination of scale lengths, we calculate a reduced statistic to indicate how well the predicted fractions 
of high- and low-a stars as a function of R and \Z\ reproduce what we see in our SEGUE sample; these values are 
indicated by the blue-shaded map in Fi gure [6l In a ddition, we calculate how well the sum of the two exponential disks 
matches the total density measured bv I Juric et al. (2008); these values arc indicated by the thin yellow contours in 
Figure ini This second constraint is not strictly a statistic, as we are comparing two smooth models. Wc normalize 
the second constraint such that it has the same 10th and 90th percentile levels as the values from the first constraint. 

We determine our best-fit scale lengths by calculating a probability for each combination of scale lengths, where the 
probability is proportional to 1"^ . All probabilities are normalized such that the total probability in the parameter 
space 1 < Lthin < 10 kpc, < ithick < 8 kpc is equal to one. The thick dashed orange contour in Figure [6] shows the 
contour that encompasses 68% of the volume under the surface defined by the probabilities. The best fit value of each 
scale length is obtained by marginalizing the probabilities over each axis and determining the peak in one dimension. 
The error bars indicate the 68% confidence interval in one dimension. This exercise yields our final results for the scale 



lengths: L 



thin 



3.C 



kpc, L 



thick 



= 1. 



3+2.1 
'-Q.5 



kpc. 



ADDITIONAL SOURCES OF ERROR 

In this section, we estimat e the errors i n our scale length estimates due to (1) random errors in the stellar parameters, 
(2) systematic errors in the I Juric et "all (|2008| ) scale lengths, and (3) c ontamina tion by halo stars. First, we estimate 
the random errors on the scale len gths usin g the same method as in iPaper ] . where we generate 500 Monte Carlo 
realizations of our data (see §6.3 of iPapeF^ for more details). In each realization, we perturb the stellar parameters 
Teff, [Fe/H], and [a/Fe] by the typical errors (200 K, 0.3 dcx and 0.1 dex, respectively). Wc find that errors in the 
stellar parameters only change the scale lengths by 0.1 kpc. 

Secondly, to estimate the errors in the assumed total density (i.e., the thin yellow contours in Figure [6]), we repeat 
the calculation, varying the I Juric et"an (j2008f ) scale lengths and scale heights by their repor ted errors (20%). The 
purpose of this exercise is to simulate the effect of systematic errors between the distances of I Juric et al.l ()2008() and 
the present work, which will cause the structural parameters to increase or decrease together. When we increase the 
lJuric et all (I2008D values by 20%, wc obtain ithin = 8.lt2'7 kpc, ithick — l-Sto'e ^pc; for a 20% decrease, we obtain 

Lthin = S.StJJ kpc, Lthick l-7tj;5 kpc. 

Lastly, to test the robustness of our results, we repeat the calculation after removing stars that may belong to the 
halo. Halo stars also have a short scale length and are enhanced in [a/Fe]. We adopt the three criteria to identify 
probable halo stars: (1) a metallicity cut that removes all stars with [Fe/H] < —0.7, (2) a kinematic cut, which removes 
all stars with < 150 km s~^, and (3) a kinematic cut that removes all stars with Vcai < 100 km s~^ to remove stars 
with the largest velocity offset relative to the projection of the local standard of rest, where Vcai = ^fl + 220 • cos6 • sinZ 
and Vr is the line-of-sight velocity measured from the SEGUE spectra. We only remove stars with Vcai < 100 km s^^ 
along lines of sight with 50 < I < 130°. We do not include the lines of sight directed toward the Galactic anticenter 
because the local standard of rest is tangent to those directions, and the projection does not give a meaningful velocity. 
For all three criteria we obtain the same scale lengths, which suggests that halo contamination docs not affect our 
scale length measurements. 
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The above analysis shows that the best-fit thick disk scale length is not significantly affected by errors in the stellar 
parameters, our assumptions of the total stellar density, and possible contamination from halo stars. The thin disk 
scale length, however, is not well-constrained because we are li mited by the lack of cov erage in R and \Z\, particularly 
in the plane of the Galaxy. Future surveys such as APOGEE f|Eisenstein et al.ll2011| ) will be able to provide stricter 
constraints on both scale lengths. 
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